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Abstract 

We discuss nuclear effects in the Paschos- Wolfenstein relationship in the context of extraction of 
the weak mixing angle. We point out that the neutron excess correction to the Paschos- Wolfenstein 
relationship for a neutron-rich target is negative and large on the scale of experimental errors 
of a recent NuTeV measurement. We found a larger neutron excess correction to the Paschos- 
Wolfenstein relationship for total cross sections than that discussed by the NuTeV collaboration. 
Phenomenological applications of this observation are discussed in the context of the NuTeV devi- 
ation. Uncertainties in the neutron excess correction are estimated. Effects due to Fermi motion, 
nuclear binding, and nuclear shadowing are also discussed in the context of total cross sections. 
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The NuTeV collaboration recently reported the results of the measurement of the weak 
mixing angle in deep inelastic neutrino and antineutrino scattering from a heavy target 
Q. The NuTeV value of sin 2 ^ = 0.2277 ± 0.0013(stat.) ± 0.0009(syst.) turned out to be 
significantly larger than that derived from a global standard model fit to other electroweak 
measurements, 0.2227 ± 0.0004 0]. The discussion of possible uncertainties and physics 
behind this discrepancy can be found in ref. [3] . However, before speculating on possible new 
physics, one should first worry about "standard" effects and uncertainties. 

A useful tool, employed by the NuTeV collaboration to derive sin 2 ^, is the Paschos- 
Wolfenstein relationship (PW) 
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where a^ c and a v CG are the deep inelastic neutrino cross sections for neutral current (NC) 
and charged- current (CC) interactions, and <r^ c and a v CQ are corresponding antineutrino 
cross sections 

In the world without heavy quarks and with exact isospin symmetry, Eq.(J^) is an exact 
relation for the isospin zero target for both the total and differential cross sections. The 
validity of this relationship is solely based on the isospin symmetry 0]. Therefore, Eq.(JIJ also 
holds for a nuclear target, provided that the nucleus is in the isoscalar state. In particular, 
this means that various strong interaction effects, including nuclear effects, must cancel out 
in the ratio (0). However, in the real world Eq.(£Q) must be corrected for the s- and c-quark 
effects. Furthermore, the targets used in neutrino experiments are usually non-isoscalar 
nuclei, and Eq.Q must also be corrected for non-isoscalarity effects, as well as for other 
nuclear effects such as Fermi motion, binding, and nuclear shadowing. 

Nuclear effects in the context of sin 2 6 W were recently discussed in refs. 0, H, H[ . It was 
suggested in ref.Q that a difference between the nuclear shadowing in NC and CC inter- 
actions may account for the discrepancy in NuTeV's measurement of sin 2 #w, although no 
specific calculation of this effect was given. The impact of nuclear effects on the extraction 
of sin 2 9w from data was discussed in ref. [£| in terms of a phenomenological parameteri- 
zation of nuclear effects in parton distributions derived from charged lepton deep inelastic 
scattering. Phenomenological studies of the difference between nuclear effects in the u- and 
(i-quark distributions in the context of the PW relationship were presented in ref. |9j. 

In this paper we address nuclear effects in the PW ratio for total cross sections. We start 
from the discussion of the non-isoscalarity correction to the PW ratio, since there appears 
to be some confusion about the sign and the magnitude of this effect in the literature. Then 
we discuss Fermi motion, nuclear binding, and nuclear shadowing corrections to the PW 
ratio. 



I. THE PW RELATIONSHIP FOR GENERIC TARGET 

We discuss (anti)neutrino deep inelastic scattering in the leading twist QCD approxima- 
tion assuming that the four-momentum transfer Q is large enough. In this approximation 
the NC and CC structure functions are given by well-known expressions in terms of quark 
and antiquark distribution functions [T3|. In order to simplify discussion of isospin effects, 
we consider the isoscalar, qo(x) = u(x) +d(x), and isovector, q\(x) = u(x) —d(x), quark dis- 
tributions (for simplicity, we suppress the explicit notation for the Q 2 dependence of parton 
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distributions). The calculation of the NC and CC cross sections, and the PW ratio in the 
QCD parton model is straightforward. A QCD radiative correction to the PW ratio for the 
total cross sections was calculated in ref.|3j. The result can be written as 
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where s 2 ^ = sin 2 6w |20|] , a s is the strong coupling, and x~ = J dxx(a a — q a ), with q a 
and q a the distribution functions of quarks and antiquarks of type a j2l|. The subscripts 
and 1 refer to the isoscalar q and isovector qi quark distributions, respectively. In the 
derivation of Eq.(J2J) we expanded in x^/xq and x~Jxq and retained only linear corrections. 
Eq.(J2J) applies to any (not necessarily isoscalar) nuclear target. We observe that the PW 
relationship (|1J) is corrected by the C-odd parts of the isovector component in the target 
(xi) and the strange and charm components of the target's sea. 

The isovector quark distribution qi vanishes in an isoscalar target, provided that the 
isospin symmetry is exact. However, a correction due to the quark-antiquark asymmetry 
in the nucleon (or nuclear) strange sea is possible. In particular, a positive x~ would move 



towards a standard model value |3j . However, available phenomenological estimates are 
controversial even in the sign of the effect: a shift in s 2 ^, estimated in ref. 3J, is —0.0026, 
however the NuTeV collaboration gives a positive shift of 0.0020 ± 0.0009 |6|. 



II. NON-ISOSCALARITY CORRECTION IN A HEAVY TARGET 

Heavy nuclei, such as iron, have unequal number of neutrons (N) and protons (Z). The 
ground state of such nuclei is not isospin zero state but rather a mixture of different isospins. 
Therefore, the isovector quark distribution is finite in such nuclei that causes a finite non- 
isoscalarity correction to structure functions and cross sections. This correction turns out 
to be large in the context of the NuTeV experiment, and, therefore, must carefully be taken 
into account. 

In order to understand this effect, we first neglect other nuclear effects and view the 
neutrino scattering off a nucleus as incoherent scattering off bound protons and neutrons at 
rest. In this approximation, the nuclear distribution of partons of type a is the sum of those 
for bound protons and neutrons 

q a/A = Zq a/p + Nq a/n = ^ (q a/p + q a/n ) + - — (q a/p - q a/n ) , (3) 

where A = Z + N. We now apply this relation to the isoscalar and the isovector quark 
distributions. Assuming exact isospin symmetry, we have q / p = q / n and qy p = —qu n . For 
nuclear parton distributions per one nucleon we then obtain 

A^qo/A = qo/p, (4a) 
A^qi/A = SNq 1/p} (4b) 

where we introduced a fractional excess of neutrons SN = (N — Z)/A. Similar equations 
can readily be written for antiquark distributions. 

Eqs.p)) suggest a negative neutron excess correction to R~ in a neutron-rich target. 
Indeed, using Eqs. P]) . we find from Eq.(J2J) that the correction is 
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where x 1 and taken for the proton. 

The NuTeV collaboration takes into account the non-isoscalarity correction in the anal- 
ysis and discusses this effect on R~ 0. However, Eq.(|5|) is different from the corre- 
sponding NuTeV equation (Eq.(7) in ref.p and Eq.(9) in ref.p). First, it must be noted 
that the NuTeV equation has the wrong sign of the SN term. Nevertheless, the NuTeV 
non-isoscalarity correction is eventually negative [11]. Therefore, in the following discussion 
we assume a correct sign in Eq.(7) in ref.p. Second, Eq.O) involves only C-odd terms, 
hence the factor x7 /xq. The corresponding factor in refs.p, |6( is (x u — Xd)/(x u + Xd), where 
x a = J dxxq a (x) for the quark distribution q a in the proton. This factor has a mixed C- 
parity and for this reason is not allowed in Eq.fjSJ). Furthermore, Eq.fjSJ) includes the a s 
correction, which was not taken into account in refs. 0| . 

In order to understand the magnitude of the non-isoscalarity correction to R~ and com- 
pare it to that in refs. 0,0, we first neglect the a s correction and compute x7 1x7 = 0.43 and 



(x u — Xd)/(x u + Xd) = 0.34 using CTEQ5 parton distributions in MS scheme |12( evaluated 
at Q 2 = 20GeV 2 , an average Q 2 in the NuTeV experiment. For the neutron excess we use 
the value 5N = 0.0567 reported by NuTeV Using s^r ~ 0.22 and collecting all factors 
in Eq.(J5J), we obtain the neutron excess correction to R~~ about —0.012. The corresponding 
correction computed using Eq.(7) in ref. is —0.0094. This correction is smaller in magni- 
tude because the factor {x u — Xd)/(x u + Xd) is smaller than the factor x7 jx7 by about 25%. 
If we simply apply the difference, —0.0026, to the NuTeV central value of sin 2 9w, we get 
a new value, 0.2251, which is now about 1.5er away from the standard model value. The 
NuTeV collaboration also reported the neutron excess correction to be —0.0080 0, [nj. m 
this case the additional shift in the NuTeV sin 2 6 W is —0.004, and the discrepancy between 
the NuTeV and the standard model value of sin 2 Qw practically disappears. However, it must 
be noted that although the non-isoscalarity correction of —0.008 was discussed in ref.|5j in 
the context of R~, it does not directly apply to the ratio of cross sections. The NuTeV 
collaboration explains that this is the shif t in the NuTeV sin 2 9w that takes into account 
experimental background and the cuts 
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The a s correction to R is also negative, but small compared to the leading term. Using 
the NLO a s = 0.2161 at Q 2 = 20GeV 2 , we find the a s correction to Br about -0.0002. 

Let us now discuss possible uncertainties in the non-isoscalarity correction. According to 
NuTeV 0], the uncertainty in this correction to sin 2 6w is 0.00005. Apparently, this is due to 
the uncertainty in the neutron excess in the target |22j. However, this is not the only source 
of uncertainties in the 5N correction. Indeed, the ratio x7/x7 is subject to theoretical 
and experimental uncertainties. This ratio depends on the set of parton distributions as 
well as on the order of perturbation theory to which the analysis is performed j23[. For 
example, the ratio x± /xq calculated with the MRST parton distributions [l3| at the same 
value of Q 2 is very similar to that of CTEQ5. However, x^/xq = 0.445 for Alekhin's 
parton distributions 0]. Furthermore, the parton distributions are known with some error. 
This error was estimated in a recent analysis in ref. 0|. Using the parton distributions 
of ref.^iJl, we have the estimation on the uncertainty in the ratio 5(x7/xq ) ~ 0.04. This 
results in the uncertainty in R~ about 0.001, which is the order of magnitude larger than the 
uncertainty in the neutron excess in the target. As an illustration, we treat this uncertainty 
as an independent theoretical uncertainty in sin 2 6w and add it in quadrature to the NuTeV 
error. Then we find some 1.2cr distance between the standard model value and the corrected 
value of NuTeV sin 2 9w discussed above. We also comment that the variations of the parton 
distributions with Q 2 introduce additional uncertainty, which is hard to access in the present 
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analysis. 



III. FERMI MOTION AND NUCLEAR BINDING CORRECTIONS 

We briefly discuss other nuclear effects on R~ . In order to sort out different effects, in 
the present discussion we do not consider explicit violations of the isospin symmetry. We 
recall that the isospin symmetry requires nuclear effects to cancel out in the R~ ratio for 
the isoscalar target. For a non-isoscalar nucleus, Eq.(j5J) was derived neglecting the effects of 
Fermi motion and nuclear binding. However, it is rather obvious that these effects and the 
non-isoscalarity correction must be considered in a unified approach. In such an approach 
Eq.® must be replaced by a convolution of quark distributions with nuclear distributions 
of bound protons and neutrons Note that in the approximation in which the proton 

and neutron distributions are identical, the Fermi motion and nuclear binding effects cancel 
out in the ratio x^/xq 24]|. Therefore, a possible correction comes through the difference 



between the proton and neutron distribution functions and is likely to be small. One source 
of the correction is a nonlinear dependence of the nucleon distribution functions on the 
number of bound particles. In heavy nuclei with Z ^ N this effect results in different 
nucleon distribution functions in the isovector and isoscalar channels even if no violation 
of the isospin symmetry is admitted [25|. Because of this effect the ratio x± /xq receives a 
correction factor of 1 — §T7M, where T is the average kinetic energy of the bound nucleon 
and M the nucleon mass jig . For the iron nucleus this factor is 0.993 that gives only a 
small correction to R~. 



IV. NUCLEAR SHADOWING EFFECT 

We now discuss nuclear shadowing effect in the context of the isoscalar and isovector 
quark distributions. The nuclear shadowing corrections to the quark distributions can be 
written as 

A^qo/A = qo/p + <Wo, (6a) 
A^qi/A = -5Nq 1/p + 5 sh gi. (6b) 

Similar equations can be written for antiquark distributions. The nuclear shadowing effect 
in the isoscalar C-even (q£ = q + q Q ) and C-odd (q$ = qo — q~o) quark distributions was 
studied in ref . [TtJ , where a different magnitude of the shadowing effect in q£ and in was 
observed. In particular, it was argued that the relative shadowing effect in the C-odd quark 
distribution is enhanced compared to that in the C-even distribution (see also ^3). m or der 
to estimate the shadowing effect in the total cross sections, we apply the approach of ref. (Tj| 
and calculate the shadowing corrections to the average quark light-cone momentum, S s hXo, 
in the C-even and C-odd distributions. For the iron nucleus the result is S s ^Xq /x^ = —0.01 
and Ss^Xq /x^ = —0.002. We observe that the shadowing effect in the total cross sections 
is significantly bigger in the C-even channel, in spite of the enhancement of the relative 
nuclear shadowing effect in the C-odd quark distributions. This, paradoxical from the first 
glance, result can be understood by observing that the C-odd cross section is saturated by 
the valence region at large x, where the nuclear shadowing effect is small. On the other 
hand, a large part of the C-even cross section comes from the sea region, where the nuclear 
shadowing effect is essential. 
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Since the isovector distribution must vanish in the isoscalar nucleus, the shadowing effect 
in the isovector channel appears at least in the order 5N. This effect has not been quan- 
titatively studied yet. In our estimates we neglect the shadowing correction to x± , leaving 
this interesting question to further studies. Then, the nuclear shadowing correction reduces 
to the renormalization of the ratio x± /x^ by the factor of 1.002 for the iron nucleus. 

Combining the nuclear shadowing effect with the Fermi motion and nuclear binding cor- 
rections to R~ , we observe a partial cancellation between these effects. For the iron nucleus 
the resulting correction factor to the ratio /xq is 0.995 that gives a negligible correction 
to R~. 

V. SUMMARY 

We have studied nuclear effects in the PW relationship for total cross sections. We ob- 
served that for a neutron-rich target (such as iron, used in the NuTeV experiment) the 
neutron excess correction to the PW relationship is negative and large on the scale of ex- 
perimental errors. We found a larger neutron excess correction than that discussed by the 
NuTeV collaboration in the context of R~ . The ratio R~ is rather sensetive to the value of 
the neutron excess correction and a consistent treatment of this correction may well explain 
at least a large part of the NuTeV deviation. 

The uncertainties in nuclear non-isoscalarity correction were discussed. We found that 
the uncertainty is dominated by the uncertainty in the parton distributions leading to the 
variation in R~ about 0.001 that should be considered as an additional source of theoretical 
uncertainty in sin 2 6w- 

We discussed the effects of Fermi motion, nuclear binding and nuclear shadowing on the 
ratio R~ for the total cross sections. We point out that nuclear effects on R~ vanish for 
the isoscalar nucleus and appear in the first or higher order in SN. For this reason these 
effects are small. Furthermore, we observed a partial cancellation between Fermi motion 
and nuclear shadowing effects in the ratio R~ . 

However, this cancellation of nuclear effects in the ratio R~ for the total cross sections 
does not, of course, mean that these effects identically cancel out in the ratio of differential 
cross sections. Since the NuTeV experiment does not measure the total cross sections, some 
nuclear effects which cancel out in R~ remain in actual experimental observables. Therefore, 
in order to clarify the impact of nuclear effects on the NuTeV result, it is necessary to 
explicitly take them into account in the analysis. 

Acknowledgments 

This work was supported in part by the RFBR project no. 00-02-17432. I am grateful 
to S.I. Alekhin, A.L. Kataev, S. Liuti, W. Melnitchouk, R. Petti, V.A. Rubakov, and G.P. 
Zeller for useful discussions and to the Theory Division of CERN, where a part of this work 
was done, for a warm hospitality. 



[1] G. P. Zeller, et al, Phys. Rev. Lett. 88, 091802 (2002). 



6 



[2] LEP Collaborations, "A Combination of Preliminary Electroweak Measurements and Con- 
straints on the Standard Model", hep-ex/0112021 
[3] S. Davidson, et al, JHEP 02, 037 (2002). 
[4] E. A. Paschos and L. Wolfenstein, Phys. Rev. D 7, 91 (1973). 
[5] K. S. McFarland, et al, in NuIntOl Proceedings, http://neutrino.kek.jp/nuint01/ 
[6] G. P. Zeller, et al, Phys. Rev. D 65, 111103 (2002). 
[7] G. A. Miller and A. W. Thomas, hep-ex/0204007 

[8] S. Kovalenko, I. Schmidt, and J.-J. Yang, Phys. Lett. B 546, 68 (2002). 

[9] S. Kumano,|hepjph70209200. 
[10] Review Of Particle Physics, Phys. Rev. D 66, 010001 (2002). 
[11] CP. Zeller, personal communications. 
[12] H. L. Lai, et al, Eur. Phys. J. C 12, 375 (2000). 
[13] A. D. Martin, et al, Eur. Phys. J. C 23, 73 (2002). 
[14] S. I. Alekhin, |hep-ph/0211096[ 

[15] S. V. Akulinichev, S. A. Kulagin, and G. M. Vagradov, Phys. Lett. B 158, 485 (1985); S. A. 

Kulagin, Nucl. Phys. A 500, 653 (1989); S. A. Kulagin, G. Piller, and W. Weise, Phys. Rev. 

C 50, 1154 (1994). 
[16] S. A. Kulagin, in preparation. 

[17] S. A. Kulagin, |hep-ph/9812532| 

[18] M. L. Mangano, et al, |hep-ph/0105155 



[19] NuTeV experiment was primarily designed to measure sin 2 6\y using the Paschos- Wolfenstein 
relationship !5]. However, it must be noted that the NuTeV collaboration does not directly 
measure the cross sections and the ratio R~ , but rather measures "ratios of experimental 
candidates within kinematic criteria and compares this to a Monte Carlo simulation" |6j. 

[20] Eq.(J2J) was derived in the tree approximation. Electroweak corrections depend on the definition 
of the weak mixing angle. To be specific, the weak mixing angle is defined here in the on-shell 
scheme, s 2 w = 1 - M^/M|. 

[21] The total cross sections involve the integration of the structure functions over the full phase 
space of x and Q 2 . Therefore a s and the moments x~ of the parton distributions are taken at 
some average scale Q 2 , which has to be chosen according to specific experimental conditions. 

[22] Using the 0.05% uncertainty in 5N and Eq.l|5]). we obtain the uncertainty in sin 2 9w about 
0.0001. 

[23] The NuTeV collaboration uses its own parton distributions obtained in the leading order fit 
to their CC data. We comment in this respect that, since the NuTeV collaboration aims at 
precize determination of the weak mixing angle, the inclusion of NLO corrections to parton 
distributions appears as a necessary improvement of the analysis. 

[24] However, it must be emphasized that this cancellation applies to the total cross sections. For 
the ratio of differential cross sections (or for the quark distributions) this effect remains finite 
even in the approximation in which the distribution functions of bound protons and neutrons 
are identical. 

[25] It is interesting to note that this effect may mask the effects of isospin violation in quark 
distribution functions. It must be also commented that possible difference between nuclear 
modifications of u- and (f-quark distributions was discussed in ref.0] in the context of the 
NuTeV deviation in terms of a phenomenological approach. 



7 



